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Abstract 

In  this  paper  we  introduce  a  model  and  an  optimization  methodology  for  terrestrial  solar  thermoelectric  generators  (STEGs).  We 
describe,  discuss,  and  justify  the  necessary  constraints  on  the  STEG  geometry  that  make  the  STEG  optimization  independent  of  indi¬ 
vidual  dimensions.  A  simplified  model  shows  that  the  thermoelectric  elements  in  STEGs  can  be  scaled  in  size  without  affecting  the  overall 
performance  of  the  device,  even  when  the  properties  of  the  thermoelectric  material  and  the  solar  absorber  are  temperature-dependent. 
Consequently,  the  amount  of  thermoelectric  material  can  be  minimized  to  be  only  a  negligible  fraction  of  the  total  system  cost.  As  an 
example,  a  Bi2Te3-based  STEG  is  optimized  for  rooftop  power  generation.  Peak  efficiency  is  predicted  to  be  5%  at  the  standard  spectrum 
AM1.5G,  with  the  thermoelectric  material  cost  below  0.05  $/Wp.  Integrating  STEGs  into  solar  hot  water  systems  for  cogeneration  adds 
electricity  at  minimal  extra  cost.  In  such  cogeneration  systems  the  electric  current  can  be  adjusted  throughout  the  day  to  favor  either 
electricity  or  hot  water  production. 

©  2012  Elsevier  Ltd.  All  rights  reserved. 
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1.  Introduction 

For  over  a  century  thermoelectric  devices  have  drawn 
little  attention  as  a  potential  technology  for  terrestrial  solar 
power  conversion  due  to  low  efficiency  and/or  complicated 
and  bulky  designs,  making  the  technology  economically 
viable  only  for  niche  solar  applications  (Telkes,  1954; 
Goldsmid  et  al.,  1980;  Dent  and  Cobble,  1982;  Rockendorf 
et  al.,  1999;  Vatcharasathien  et  al.,  2005;  Li  et  al.,  2010). 
Solar  thermoelectric  generators  (STEGs)  have  also  been 
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designed  and  optimized  for  space  applications  due  to  their 
advantages  of  reliability  and  the  capability  to  withstand 
high  incident  solar  radiation  (Fuschillo  et  al.,  1966;  Scher- 
rer  et  al.,  2003)  However,  constraints  and  designs  for  earth- 
based  STEGs  are  different  and  have  not  been  successfully 
optimized  for  large-scale  and  potentially  cost-effective 
deployment  on  rooftops  or  for  the  integration  in  solar 
hot  water  systems.  The  most  recent  experimental  work  on 
STEG  technology  showed  promising  results  for  those  types 
of  applications  (Kraemer  et  al.,  2011). 

The  simplest  design  of  a  STEG  consists  of  a  unicouple 
thermoelectric  generator  (TEG)  sandwiched  between  a 
solar  absorber  that  absorbs  the  incident  solar  radiation 
and  a  heat  sink  that  removes  the  excess  heat  to  maintain 
the  temperature  difference  across  the  unicouple  TEG 
(Fig.  1).  Typically  a  unicouple  STEG  includes  optical 
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Fig.  1.  Concept  of  solar  thermoelectric  power  conversion,  (a)  A  solar-driven  TEG  cell  consists  of  a  possible  optical  concentrator,  a  solar  absorber 
(wavelength-selective  surface  on  a  copper  substrate)  which  absorbs  the  solar  radiation  and  thermally  concentrates  the  heat  onto  the  n/p-type  unicouple 
TEG.  The  opposite  side  of  the  TEG  is  attached  to  copper  electrodes  that  are  used  to  connect  the  electrical  circuit  to  extract  the  electrical  power  from  the 
system;  (b)  conceptual  block  diagram  illustrating  a  solar  thermoelectric  generator  (STEG)  with  its  model  components,  the  corresponding  input  parameter 
and  the  logic  to  calculate  the  heat  flows  and  the  electrical  power  output  of  the  STEG. 


components,  such  as  a  glass  enclosure  to  maintain  a  vac¬ 
uum  and/or  an  optical  concentrator  system  to  increase 
the  incident  solar  radiation  on  the  absorber  to  reduce  heat 
losses.  The  use  of  thermal  concentration  to  focus  absorbed 
solar  radiation  flux  via  lateral  heat  conduction  within  the 
absorber  onto  the  thermoelectric  generator  has  been  pro¬ 
posed  in  several  studies  (Telkes,  1954;  Fuschillo  et  al., 
1966;  Goldsmid  et  al.,  1980;  Rowe,  1981;  Kraemer  et  al., 
2011;  Chen,  2011).  The  thermal  concentration,  Ctb,  is 
defined  as  the  geometric  ratio  of  the  absorber  surface  area, 
Aabs,  to  the  TEG  cross-sectional  area,  Ateg.  Just  like  optical 
concentration  (Duffle  and  Beckman,  2006),  the  defined 
thermal  concentration  is  not  to  be  confused  with  the  actual 
energy  flux  concentration  which  is  also  dependent  on  the 
properties  and  temperature  of  the  solar  absorber.  Using 
thermal  concentration  reduces  the  amount  of  required 
thermoelectric  material,  and  also  increases  the  TEG  effi¬ 
ciency  by  creating  a  higher  absorber  temperature.  How¬ 
ever,  the  thermal  losses  from  the  absorber  become 
significant  with  increasing  absorber  temperature,  suggest¬ 
ing  the  existence  of  an  optimal  operational  absorber  tem¬ 
perature  which  balances  the  TEG  efficiency  gains  with 
the  radiation  losses  (Telkes,  1954;  Rowe,  1981;  Kraemer 
et  al.,  2011;  Chen,  2011).  A  conceptual  picture  of  the  opti¬ 
mization  of  a  STEG  follows.  A  unicouple  STEG  comprises 
an  optothermal  system  to  convert  radiation  into  a  heat 
flux,  and  a  unicouple  TEG  to  directly  convert  the  heat  flux 
into  electricity.  The  efficiency  of  a  given  optothermal  sys¬ 
tem  under  given  ambient  conditions  is  a  function  of  the 
temperature  of  the  absorber,  as  material  properties  and 
radiation  losses  are  a  function  of  temperature.  The  effi¬ 
ciency  of  a  unicouple  TEG  of  given  material  properties 
and  ambient  conditions  is  also  a  function  of  the  tempera¬ 
ture  of  the  absorber.  Essentially,  this  absorber  temperature 
is  proportional  to  the  product  of  the  thermal  resistance  of 
the  unicouple  TEG,  and  the  heat  flow  into  it.  The 


geometric  parameter  which  determines  the  heat  flow  into 
the  unicouple  TEG  is  the  absorber  area;  the  geometric 
parameter  which  determines  the  thermal  resistance  of  the 
unicouple  TEG  is  the  ratio  (L/Ateg)  where  L  is  the  length 
of  the  thermoelectric  elements  and  Ateg  is  the  cross-sec¬ 
tional  area  of  the  unicouple  TEG.  Therefore  the  geometric 
parameter  which  controls  the  absorber  temperature,  and 
thus  the  efficiency  of  a  STEG,  is  the  product 
CthL  =  Aabs(L/ Ateg).  This  simple,  intuitive  result  was  math¬ 
ematically  derived  by  Chen  (2011)  for  a  STEG  with  large 
thermal  concentration,  but  only  holds  as  long  as  certain 
conditions  are  met,  which  will  be  explained  in  this  paper. 
Prior  to  the  analysis  by  Chen  (2011)  a  STEG  design  with 
large  optical  concentration  was  presented  by  Rowe  (1981) 
based  on  similarly  simplified  models  assuming  constant 
thermoelectric  and  solar  absorber  material  properties, 
geometry-independent  radiative  heat  losses,  no  electrical 
contact  resistance,  and  no  thermal  resistance  in  the  solar 
absorber.  Also,  the  inevitable  variations  of  the  incident 
solar  radiation  and  cold-junction  temperature  of  the  STEG 
over  the  course  of  the  day  and  the  effect  on  the  device  opti¬ 
mization  have  not  been  addressed  by  previous  publications. 

In  this  paper  we  describe  a  STEG  model  accounting  for 
various  effects  neglected  by  previous  publications  such  as 
Rowe  (1981)  and  Chen  (2011)  and  discuss  in  detail  when 
and  how  those  effects  will  influence  the  STEG  performance 
and  the  STEG  geometric  optimization  parameter.  We  show 
that  only  within  certain  limits  the  assumptions  regarding 
contact  resistance,  thermal  resistance,  and  geometry- 
dependent  losses  are  valid.  Simplifying  the  detailed  model 
by  introducing  the  justified  simplifications  allows  to  reduce 
the  computational  time  significantly  without  losing  accu¬ 
racy.  The  model  accounts  for  the  temperature  dependence 
of  the  thermoelectric  materials  and  solar  absorber  proper¬ 
ties  while  neglecting  electrical  contact  resistances,  the  non¬ 
uniformity  of  the  absorber  temperature,  and  the  radiation 
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exchange  with  the  thermoelectric  elements.  Using  this 
model  we  will  discuss  the  influence  of  the  inevitable  varia¬ 
tion  of  the  solar  intensity  over  the  course  of  a  day  and  the 
cold-junction  temperature  on  the  STEG  performance  and 
the  STEG’s  geometric  optimization  parameter. 

We  also  use  this  model  to  analyze  the  cogeneration 
application  by  integrating  a  STEG  in  a  solar  hot  water  sys¬ 
tem  providing  electricity  and  hot  water.  Even  though  this 
paper  is  mainly  focused  on  STEGs  with  small  or  no  optical 
concentration  and  large  thermal  concentration,  the  intro¬ 
duced  optimization  methodology  and  model  can  also  be 
used  for  STEGs  with  large  optical  concentration. 

2.  General  model 

2.1.  Geometry 

The  modeled  STEG  design  (Fig.  la)  consists  of  a  TEG 
unicouple  (i.e.  a  pair  of  n-  and  p-type  doped  Bi2Te3  ele¬ 
ments)  sandwiched  between  a  wavelength-selective  solar 
absorber  (Kennedy,  2002)  and  two  polished  copper  elec¬ 
trodes,  surrounded  by  a  glass  vacuum  enclosure.  The  sys¬ 
tem  can  be  placed  in  the  focal  point  of  an  optical 
concentrator  system  to  increase  the  operational  tempera¬ 
ture  difference  and  conversion  efficiency  of  the  STEG  (Kra¬ 
emer  et  al.,  201 1).  However,  this  work  is  mainly  focused  on 
non-optical  concentrating  STEGs  with  large  thermal  con¬ 
centration.  For  the  ease  of  the  STEG  fabrication  the  legs 
are  of  equal  length  but  can  have  different  cross-sectional 
areas.  The  solar  absorber  is  a  thin  copper  substrate  coated 
with  a  wavelength-selective  surface  on  the  top  side  which 
absorbs  the  whole  wavelength  spectrum  of  the  incident 
solar  radiation.  The  copper  substrate  then  conducts  the 
heat  laterally  to  the  unicouple  TEG  which  transports  the 
concentrated  heat  to  the  cold  junction  while  converting 
part  of  the  heat  into  electricity.  A  heat  sink  removes  the 
heat  from  the  cold  junction,  where  it  is  rejected  to  the  sur¬ 
roundings  or  used  for  cogeneration. 

2.2.  Input  parameters 

The  inputs  to  the  optics  model  are  the  effective  incident 
radiation  on  the  plane  of  the  aperture,  qsof,  the  optical  effi¬ 
ciency;  and  possible  geometric  optical  concentration.  The 
thermal  model  of  the  solar  absorber  requires  the  tempera¬ 
ture-dependent  emittance  of  the  selective  surface;  the 
absorber  geometry;  and  the  thermal  conductivity  of  the 
absorber  substrate.  The  TEG  model  requires  the  geometry 
of  the  different  thermoelectric  elements;  the  temperature- 
dependent  thermoelectric  material  properties;  the  external 
electrical  load;  and  the  electrical  contact  resistance  at  the 
ends  of  the  thermoelectric  elements.  The  cold-side  model 
requires  the  thermal  conductance  of  the  cold  side.  All  four 
models  also  rely  on  the  ambient  temperature,  and  the 
geometry-related  view  factors  between  all  the  parts  of  the 
device.  For  the  radiative  heat  balance  of  the  back  side  of 
the  solar  absorber  with  its  surroundings  (cold  side, 


thermoelectric  elements,  gap  between  absorber  and  cold 
side)  the  emittance  value  of  each  participating  surface  must 
be  know. 

2.3.  Governing  equations 

The  STEG  efficiency  is  defined  as  the  ratio  between  the 
electrical  power  output,  Wei,  and  the  incident  solar  radia¬ 
tion  power,  Qsoi,  which  is  the  product  of  the  effective  inci¬ 
dent  solar  radiation  flux,  qsoh  and  the  aperture  area,  Aap,  of 
the  STEG  (Eq.  (1)) 


While  the  denominator  of  Eq.  (1)  is  straightforward,  the 
electric  power  must  be  calculated  by  using  the  submodels 
(Fig.  lb).  Without  any  simplifications,  all  four  submodels 
must  be  solved  simultaneously  because  of  the  radiation 
exchange  between  them.  We  will  describe  the  submodels 
in  detail  below. 

The  optical  model  is  a  general  model  that  lumps 
together  the  absorptance,  at,  of  the  solar  absorber,  the  glass 
transmittance,  t,  possible  mirror  reflection  losses,  the  inter¬ 
cept  factor,  and  the  incident  angle  modifier  to  an  optical 
efficiency,  qopt  =  Qabs/QSoi  (Duffie  and  Beckman,  2006). 
This  optical  efficiency  is  defined  as  the  ratio  of  the  absorbed 
solar  radiation  to  the  solar  radiation  intercepted  by  the 
total  aperture.  In  the  following  simulations  the  solar 
absorptance  is  assumed  to  be  independent  of  temperature, 
which  is  a  reasonable  assumption  for  commercially-avail- 
able  multilayer  thin-film  selective  surfaces  at  temperatures 
below  250  °C.  In  the  case  of  an  optically-concentrating 
STEG  system  the  optically  concentration,  Copt,  is  the  ratio 
of  the  aperture  area  to  the  absorber  area,  Aabs.  The 
absorbed  solar  radiation  heat  flux  can  then  be  obtained 
from  Eq.  (2). 

qabs  =  Coptqoptqsoi  (2) 

The  optical  performance  of  collectors  such  as  the  inter¬ 
cept  factor  and  the  incident  angle  modifier  are  specific  to 
the  optical  design  and  beyond  the  scope  of  this  work;  more 
detailed  discussions  can  be  found  in  the  literature  (Duffie 
and  Beckman,  2006).  For  the  following  simulation  results 
the  optical  concentration,  the  intercept  factor,  and  the  inci¬ 
dent  angle  modifier  are  set  to  1. 

The  thermal  model  begins  with  the  calculation  of  the 
temperature  distribution  within  the  solar  absorber.  A  large 
lateral  temperature  drop  within  the  absorber  is  not  favor¬ 
able  for  the  system  efficiency.  The  square-shaped  solar 
absorber  is  approximated  as  a  circular  disk  of  equivalent 
surface  area  and  is  modeled  as  a  radial  fin  with  an  absorbed 
radiation  power,  Qabs,  distributed  evenly  over  it  (Fig.  2a). 
The  closely-spaced  thermoelectric  elements  are  approxi¬ 
mated  as  one  cylindrical  TEG  with  an  equivalent  cross-sec¬ 
tional  area  at  a  uniform  hot-junction  temperature.  The 
largest  radial  temperature  drop  within  the  solar  absorber 
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Fig.  2.  STEG  submodels,  (a)  Geometric  model  for  STEG  radiation  calculations;  (b)  energy  balance  over  an  absorber  segment  for  the  radially  discretized 
solar  absorber;  (c)  2D  Schematic  of  a  STEG  showing  the  merged  unicouple  thermoelectric  generator  (TEG)  with  radiative  heat  losses;  EbiSS,  Ebl,  Eb2,  and 
Eb 3  are  emissive  powers  of  the  surfaces;  Qieg,h  and  QieStC  is  the  heat  absorbed  at  the  hot  junction  and  rejected  at  the  cold  junction;  (d)  ID  discretization  of  a 
thermoelectric  element  and  energy  balance  of  a  segment  for  iterative  numerical  calculations. 


occurs  close  to  the  unicouple  TEG  due  to  the  large  heat 
flux  and  thermal  resistance.  For  the  STEG  design  consid¬ 
ered  in  this  work  with  an  absorber  area  much  larger  than 
the  cross-sectional  area  of  the  thermoelectric  elements  it 
is  reasonable  to  assume  radial  heat  flow  close  to  the  uni¬ 
couple  TEG.  Therefore  the  assumption  of  a  circular  solar 
absorber  of  equivalent  area  to  estimate  its  temperature  dis¬ 
tribution  is  reasonable.  The  top  side  of  the  absorber  has  a 
temperature-dependent  emittance,  ess.  The  back  side  of  the 
absorber  experiences  radiative  heat  transfer  with  the  cold 
side,  the  thermoelectric  elements,  and  the  environment 
through  the  gap  between  the  absorber  and  the  cold  side. 
Within  the  absorber,  the  heat  is  conducted  radially  by  the 
copper  substrate  of  thickness  t  and  thermal  conductivity 
kabs ■  The  absorber  is  discretized  radially.  From  the  energy 
balance  over  one  elemental  ring  of  the  absorber  (Fig.  2b) 
Eq.  (3)  is  obtained  which  then  can  be  solved  for  the  temper¬ 
ature  non-uniformity  within  the  solar  absorber  (Fig.  2a) 
using  a  finite  difference  scheme. 


dTabs 

dr 


d2Tabs 
T  dr2 


(2  r  +  dr) 
2kabst 


[qabs-ess(T(r)) 


x  ^(^1  -  Tl)  -  qradJk] 


(3) 


Too;  the  third  term  represents  the  net  radiative  heat  flux 
from  the  back  side  of  the  solar  absorber,  qrad,b-  The  net  radi¬ 
ative  heat  flux  couples  the  thermal  absorber  model,  the 
TEG  model,  and  the  cold-side  model  (Fig.  lb  and  2c), 
and  is  described  with  the  thermoelectric  model  below.  Eq. 
(3)  is  solved  numerically  with  the  hot-junction  temperature 
as  the  boundary  condition  at  the  base  of  the  fin, 
T(r  =  rteg)  =  Tbase  —  Ttegji  and  with  the  absorber  edge  radi¬ 
ation  flux,  qrad.e  =  ecuo(T*bs(r  =  Rabs)  —  T4x),  as  a  Neumann 
boundary  condition  for  the  absorber  with  radius,  Rabs.  The 
emittance  of  the  polished  copper  substrate,  ec„,  is  assumed 
to  be  0.03  (Mills,  1999).  The  absorber  is  assumed  to  be  at 
a  uniform  hot-junction  temperature  in  the  region  (r  <  rteg) 
where  the  absorber  is  thermally  attached  to  the  TEG.  This 
is  a  necessary  assumption  for  the  later-introduced  ID  dis¬ 
cretization  model  of  the  thermoelectric  elements.  A  detailed 
discussion  in  Section  3.3  will  justify  this  assumption  for 
STEG  designs  optimized  according  to  the  methodology 
proposed  by  this  paper.  In  other  STEG  designs  this  assump¬ 
tion  might  not  be  valid  and  will  lead  to  an  overestimation  of 
the  STEG  performance.  The  temperature  distribution  of  the 
absorber  is  determined  and  the  heat  flux  at  the  hot-junction 
of  the  TEG  can  be  obtained  from  Eq.  (4). 


The  first  term  in  the  brackets  represents  the  absorbed  \  a 

incident  solar  radiation;  the  second  term  accounts  for  the  dtegfi  =  \<labs  ~  ^  A,g„ G.s ( TabAr) )  ( Tahs (r)  —  T^jdA 
net  radiative  heat  flux  from  the  selective  surface  of  the  solar 

absorber  at  a  temperature  Tabs{R)  to  the  glass  enclosure  _  Qrad’h  +  1  ^ 

treated  as  a  blackbody  surface  at  ambient  temperature  J 


(4) 
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The  heat  flux  into  the  TEG  becomes  the  boundary  con¬ 
dition  for  the  TEG  model.  The  first  term  in  the  brackets  is 
the  absorbed  radiation  heat  flux,  the  second  term  is  the 
integrated  net  radiative  heat  loss  from  the  wavelength- 
selective  surface  to  the  blackbody  surroundings,  and  the 
third  term  consists  of  the  net  thermal  radiation  heat  flow 
from  the  copper  back  side,  Qrad.b  =  (Aabs  -  Ateg)qradh,  and 
from  the  edges,  Qrad,e  =  {2nRabst)qra<i,e,  of  the  solar  absor¬ 
ber.  The  edge  heat  losses  can  be  minimized  by  choosing  a 
thin  copper  substrate.  However,  reducing  the  thickness  of 
the  copper  substrate  increases  the  radial  thermal  resistance 
and  will  result  in  larger  radial  temperature  gradients  within 
the  absorber  as  will  be  discussed  in  Section  3.3. 

The  TEG  model  is  a  discretized  model  of  the  thermo¬ 
electric  elements  which  incorporates  thermoelectric  effects 
as  well  as  contact  resistance,  radiation  exchange  of  the  legs 
with  the  surroundings,  and  the  external  electrical  load.  The 
model,  shown  in  Fig.  2d,  is  described  in  Eq.  (5). 


1  dE 
A  ieg  dx 


^legdradix) 

A  leg 

d(iST)  1  d(fii) 

dx  e  dx 


(5) 


This  differential  equation  is  a  one-dimensional  energy 
flux  balance  over  a  segment  of  a  discretized  thermoelectric 
element  at  position  x,  where  the  left  hand  side  represents 
the  change  in  total  energy  dE  in  the  segment  of  length  dx 
divided  by  the  cross-sectional  area,  A!eg,  of  the  thermoelec¬ 
tric  element.  This  change  in  energy  dE/dx  is  equal  to  the  net 
radiative  heat  flux,  qrad,  to  the  thermoelectric  element  seg¬ 
ment  at  position  x  from  its  surroundings  multiplied  by 
the  leg  equivalent  perimeter,  3Pieg,  of  the  thermoelectric  ele¬ 
ment.  A  detailed  multi-body  radiative  heat  transfer  model 
with  the  emissive  powers  Eb^s,  Eh  u  Eb2,  and  Eh3  of  each 
contributing  surface  (Fig.  2c)  determines  the  radiative  heat 
loss  from  the  unicouple  STEG,  including  the  radiation 
from  the  edges  of  the  solar  absorber  and  from  the  thermo¬ 
electric  elements.  In  the  model  the  two  thermoelectric  ele¬ 
ments  are  merged  to  one  cylindrical  thermoelectric 
element  with  an  equivalent  surface  area 
ASieq  =  3 L(y/Ap  +  \[A^n)  to  simplify  the  view  factor  calcula¬ 
tions  (Modest,  2003)  where  Ap  and  An  are  the  cross-sec¬ 
tional  areas  of  the  p-  and  n-type  thermoelectric  elements. 
The  equivalent  surface  area  is  based  on  the  assumption  that 
the  net  radiative  heat  transfer  from  the  side  walls  of  the 
thermoelectric  elements  facing  each  other  is  negligible 
because  the  closely  spaced  n-  and  p-type  thermoelectric  ele¬ 
ments  have  similar  temperature  profiles.  Consequently,  the 
facing  side  walls  mainly  radiate  to  one  another  at  similar 
temperatures,  which  reduces  the  net  radiative  heat  transfer 
significantly.  Therefore  the  facing  side  walls  act  as  radia¬ 
tion  shields  for  one  another.  The  corresponding  effective 
leg  perimeter  is  & ieg  =  3  y/A^g.  The  thermoelectric  elements 
are  closely  spaced  in  order  to  obtain  mechanically-reliable 
unicouple  STEGs  with  stable  electrical  and  thermal 


contacts.  The  thermoelectric  elements  are  soldered  to  the 
solar  absorber  which  expands  and  contracts  due  to  thermal 
cycling  during  operation.  Consequently,  the  thermal 
expansion/contraction  will  apply  shear  stresses  on  the  sol¬ 
der  joints.  The  larger  the  distance  between  the  thermoelec¬ 
tric  elements  the  larger  would  be  the  displacement  and  thus 
the  shear  stresses  that  cause  the  degradation  of  the  contacts 
and  ultimately  the  failure  of  the  joints. 

The  thermoelectric  elements  are  discretized  in  small  iso¬ 
thermal  sections  and  the  view  factors  can  be  determined  with 
view  factor  algebra  (Modest,  2003)  by  approximating  the 
square-shaped  absorber  and  cold  side  as  two  finite  disks  at 
uniform  temperature  of  equivalent  surface  area  with  radius 
Rabs  (Fig.  2a).  In  order  to  take  into  account  the  increase  in 
radiation  losses  from  the  solar  absorber  back  side  due  to 
the  possible  non-uniformity  of  its  temperature,  the  tempera¬ 
ture  is  averaged  over  the  equivalent  radius  of  the  absorber. 
After  all  view  factors  are  determined  the  net  radiative  heat 
flux  from  each  surface  can  be  obtained  with  the  net  radiation 
method  (Modest,  2003).  The  emittance  of  the  back  side  and 
edges  of  the  solar  absorber  (polished  copper),  and  the  cold 
side  (polished  copper)  is  assumed  to  be  temperature  indepen¬ 
dent  with  ecu  =  0.03  (Mills,  1999).  The  emittance  of  the  ther¬ 
moelectric  elements  is  assumed  to  be  constant  with  ete  =  0.5. 
For  more  accurate  performance  predictions  of  a  unicouple 
STEG  with  long  thermoelectric  elements  the  temperature 
dependence  of  the  emittance  should  be  included.  However, 
for  a  high-performance  and  cost-effective  STEG  design  with 
large  thermal  concentration  and  short  thermoelectric  ele¬ 
ments  the  accurate  knowledge  of  their  radiative  properties 
is  not  important  as  the  analysis  will  show.  The  first  term  on 
the  right  hand  side  of  Eq.  (5)  represents  the  change  in  heat 
flux  related  to  the  entropy  carried  by  the  electrical  current 
flux,  The  second  term  is  the  change  in  energy  flux  due  to 
the  change  in  the  electrochemical  potential,  Vju,  of  the  energy 
carriers  divided  by  the  electron  charge  (-e): 
1  /(— e)  Wji  =  —SWT  —  pi,  where  S  is  the  temperature-depen¬ 
dent  Seebeck  coefficient  of  the  material,  and  p  is  the  electrical 
resistivity.  This  equality  can  be  derived  from  the  theory  of 
irreversible  thermodynamics  and  is  referred  to  as  the  general¬ 
ized  Ohm’s  law  (Domenicali,  1953).  The  third  term  is  the 
change  in  the  conducted  heat  with  the  thermal  conductivity 
k(x).  An  equation  for  the  local  thermoelectric  heat  flux  can 
be  obtained  from  the  entropy  flux  Js(x)  =  iS(x )  -  ^4  ^  at 
position  x  inside  the  thermoelectric  element.  The  Seebeck 
coefficient  S(x)  is  defined  as  the  entropy  carried  per  electron 
charge  (Domenicali,  1953).  Multiplying  the  entropy  flux  by 
the  local  temperature  yields  the  local  thermoelectric  heat 
flux. 

qleg(x)=i(ST)x-(k^j  (6) 

The  differential  equations  Eqs.  (5)  and  (6)  can  be 
arranged  into  a  coupled  set  of  first-order  equations 
(Mahan,  1991;  Hogan  and  Shih,  2006)  and  are  shown 
below,  including  the  net  thermal  radiative  heat  transfer  of 
the  thermoelectric  element  side  walls. 


D.  Kraemer  et  al.  I  Solar  Energy  86  (2012)  1338-1350 


1343 


dT  _i(ST)x-qleg(x) 
dx 

3  qradjx) 


(7) 

(8) 


Those  coupled  first-order  equations  are  solved  itera¬ 
tively  (Hogan  and  Shih,  2006;  Buist,  1995)  for  the  n-  and 
p-type  elements  (with  known  temperature-dependent  prop¬ 
erties)  for  the  optimal  thermoelectric  heat  flux  and  temper¬ 
ature  distribution  corresponding  to  the  optimum  current 
densities  for  a  specified  TEG  geometry  (lengths  and 
cross-sectional  areas  of  thermoelectric  elements).  The  ther¬ 
moelectric  elements  are  chosen  to  have  the  same  length. 
Due  to  the  difference  in  the  n-  and  p-type  thermoelectric 
material  properties,  the  cross-sectional  area  of  the  thermo¬ 
electric  elements  are  different  and  the  cross-sectional  area 
ratio  must  be  optimized  during  the  iteration  process.  The 
boundary  conditions  for  the  differential  equations  of  both 
thermoelectric  elements  are  the  same  fixed  cold-  and  hot- 
junction  temperatures  that  uniquely  determine  the  optimal 
current,  the  temperature  and  heat  flux  distribution  for  the 
specified  unicouple  TEG  geometry.  Electrical  contact  resis¬ 
tances  are  included  in  the  calculations  by  adding  a  small 
segment  with  corresponding  electrical  resistivity  and  zero 
Seebeck  coefficient  to  the  hot  and  cold  end  of  the  thermo¬ 
electric  elements.  From  the  determined  energy  flows  enter¬ 
ing  and  leaving  the  TEG,  the  electrical  power  output  of  the 
STEG  is  found  with  the  following  equation. 


Wei  =  A legqleg,h  -  {Qteg,n  +  QtegP}X=t  + 

x  /  q radix)  dx 
Jo 


(9) 


The  first  term  on  the  right-hand  side  of  Eq.  (9)  is  the 
TEG  heat  flux  at  the  hot  junction,  qtegih,  which  is  the 
heat  flux  absorbed  and  concentrated  by  the  solar  absor¬ 
ber.  The  term  in  braces  is  the  heat  rejected,  Qleg,c  at  the 
cold  junction  (x  =  L)  by  the  n-  and  p-type  thermoelectric 
elements.  The  last  term  is  the  net  radiative  heat  transfer 
of  the  TEG  with  its  surroundings  with 
Pteg  —  3(,V^p  +  V^Q  as  the  perimeter  and  qmd(x)  as  the 
net  radiative  heat  flux  of  the  unicouple  TEG.  The  excess 
heat,  Qexcess,  that  must  be  removed  on  the  cold  side  to 
maintain  the  optimal  temperature  difference  across  the 
TEG  can  then  be  defined  with  Eq.  (10). 

Qexcess  =  Qteg,c  +  Qrad,  =  U[TeJ  ~  Tc)  (10) 

The  cold  side  is  modeled  as  a  heat  sink  to  the  ambient 
or  to  a  working  fluid  with  an  overall  heat  transfer  coef¬ 
ficient  (HTC),  U,  which  is  assumed  to  be  infinitely  large 
for  some  simulations  to  simplify  the  calculations  by  using 
a  cold-junction  temperature,  Tcj,  which  is  equal  to  the 
fixed  cold-side  temperature,  Tc.  The  net  radiative  heat 
transfer,  Qrad,c,  between  the  cold  side  and  the  other  parts 
of  the  system  is  calculated  from  the  aforementioned  radi¬ 
ation  model. 


2.4.  Optimization  for  fixed  incident  radiation  flux 

With  the  introduced  models  the  geometry  of  the  unicou¬ 
ple  STEG  (Fig.  la)  must  be  optimized  for  optimal  hot- 
junction  temperature  which  yields  highest  electric  power 
output  (Eq.  (9))  and  corresponding  highest  STEG  effi¬ 
ciency  as  discussed  in  the  introduction.  The  geometrical 
optimization  is  performed  for  a  constant  incident  solar 
radiation  flux  and  a  fixed  cold-side  temperature.  Further¬ 
more,  some  constraints  are  enforced  on  the  optimization. 
The  length  and  the  current  of  both  thermoelectric  elements 
are  the  same.  Additionally,  the  cold-junction  and  hot-junc- 
tion  temperatures  of  both  thermoelectric  legs  are  set  to  be 
equal.  The  temperature  and  the  heat  flow  at  the  hot  junc¬ 
tion  are  introduced  as  matching  boundary  conditions  for 
the  solar  absorber  and  the  hot  junction  of  the  unicouple 
TEG.  Under  those  constraints  all  dimensions  of  the  uni¬ 
couple  STEG  (length  and  cross-sectional  areas  of  thermo¬ 
electric  elements,  and  the  absorber  size)  are  varied,  and 
load  matching  for  each  set  of  geometries  is  performed  to 
find  the  optimal  geometry. 

3.  Discussion 

The  performance  of  a  STEG  is  negatively  affected  by 
contact  resistance  losses,  the  temperature  non-uniformity 
of  the  absorber,  the  radiation  losses  out  the  gap  between 
the  absorber  and  cold  side,  and  the  radiation  losses  from 
the  sidewalls  of  the  TEG.  These  four  loss  mechanisms  are 
all  geometry-dependent,  and  are  included  in  the  introduced 
detailed  model.  However,  those  losses  can  be  minimized  if 
the  system  geometry  meets  certain  constraints.  In  the  first 
part  of  this  section,  the  necessary  constraints  for  each  of 
these  loss  mechanisms  are  developed.  When  the  STEG 
design  meets  these  constraints,  the  STEG  can  be  modeled 
in  a  much  simpler  manner.  In  the  second  part  of  this  sec¬ 
tion  a  STEG  optimization  methodology  based  on  this  sim¬ 
plified  model  is  introduced  and  discussed.  For  all 
simulations  the  STEG  comprises  an  evacuated  glass  tube 
with  an  effective  transmittance  of  90.5%;  a  commercial 
wavelength-selective  solar  absorber  with  known  tempera¬ 
ture-dependent  properties,  a  copper  substrate  thermal  con¬ 
ductivity  of  380  W/mK,  and  a  thickness  of  0.2  mm; 
nanostructured  Bi2Te3  compound  thermoelectric  materials, 
as  described  by  Poudel  et  al.  (2008);  AM1.5G  (ASTM, 
2008)  solar  flux  striking  the  glass  tube;  and  an  ambient 
and  cold-junction  temperature  of  25  °C  unless  stated 
differently. 

3.1.  Effect  of  electrical  contact  resistance 

In  order  to  obtain  the  best  performance  of  a  STEG,  the 
relative  electrical  contact  resistance  must  be  minimized.  To 
gauge  how  the  geometry  affects  the  contact  resistance 
losses,  the  other  geometry-dependent  losses  are  turned  off 
in  the  simulation.  This  is  accomplished  in  three  steps:  giv¬ 
ing  the  absorber  infinite  thermal  conductivity  eliminates 
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temperature  gradients  in  the  absorber;  setting  to  zero  the 
emittance  of  the  thermoelectric  elements  eliminates  radia¬ 
tion  losses  from  the  thermoelectric  elements;  and  setting 
to  1  the  view  factor  from  the  absorber  to  the  cold  side 
(the  infinite  parallel  plate  assumption)  prevents  the  absor¬ 
ber  back  side  losses  from  being  geometry-dependent. 
Fig.  3a  shows  the  effect  of  different  electrical  contact  resis¬ 
tances  on  the  efficiency  of  a  unicouple  STEG  under  these 
conditions.  These  assumptions  yield  an  upper  limit  for 
the  STEG  efficiency  of  5.1%.  The  results  show  that  the 
reduction  of  the  STEG  efficiency  due  to  electrical  contact 
resistance  becomes  larger  with  decreasing  thermoelectric 
element  length,  which  has  been  pointed  out  before  for 
other  thermoelectric  devices  (Goldsmid,  1986).  The  mini¬ 
mum  electrical  contact  resistance  for  a  bismuth  telluride 
to  metal  interface  is  estimated  to  be  on  the  order  of 
10-11  Qcm2  (Silva  et  al.,  2004).  Gupta  et  al.  (2010)  exper¬ 
imentally  demonstrated  electrical  contact  resistances  below 


Fig.  3.  The  effect  of  the  system  losses  of  a  simple  flat-plate  unicouple 
STEG  geometry  on  the  STEG  performance,  (a)  Effect  of  the  electrical 
contact  resistance  rc:  1 0  7  £2  cm2  solid  line,  5  x  1 0  7  ft  cm2  dashed  line, 
KT6£2cm2  dash-dotted  line;  (b)  effect  of  geometry-dependent  radiative 
heat  losses  for  a  unicouple  STEG  with  constant  finite  parallel  plates  and 
absorber  size  of  91 1  mm2  including  (solid  line)  and  excluding  (dashed  line) 
thermoelectric  element  radiation  losses  compared  to  a  unicouple  STEG 
assuming  infinite  parallel  plates  excluding  the  thermoelectric  elements 
from  the  heat  balance  (dash-dotted  line);  the  cross-sectional  area  of  the 
unicouple  TEG  is  altered  to  find  optimum  Ct/,L;  and  (c)  effect  of  the  non- 
uniformity  of  solar  absorber  temperature;  for  each  length  of  the 
thermoelectric  element  the  optimal  CthL  is  found  in  one  case  by  varying 
the  cross-sectional  area  of  the  TEG  ( A,eg )  while  keeping  the  absorber  area 
(■Aabs)  constant  (911  mm2)  and  in  the  other  case  by  changing  Aabs  while 
keeping  Ateg  constant  (3.645  mm2). 


10-7  Q  cm2  for  Bi2Te3  sputter-coated  with  Ni  or  Co.  If  the 
electrical  contact  resistance  is  as  small  as  10-7  Q  cm2  the 
effect  on  the  STEG  performance  is  negligible  even  for  ther¬ 
moelectric  elements  as  short  as  0.1  mm  (Fig.  3a).  Other 
plating  processes  can  result  in  higher  electrical  Fig.  3contact 
resistances  which  then  cannot  be  neglected  in  the  design 
optimization  and  will  limit  the  minimum  value  of  the  ther¬ 
moelectric  element  length. 

3.2.  Geometry-dependent  radiation  losses 

In  a  STEG  cell  design  (Fig.  lb),  the  STEG  geometry 
affects  the  radiation  exchange  between  the  TEG,  the  cold 
side,  the  back  side  of  the  absorber,  and  the  environment. 
Ideally  the  sidewalls  of  the  TEG  would  not  contribute  to 
the  radiative  heat  balance,  and  the  absorber  and  cold  side 
would  be  infinite  parallel  plates  where  no  radiation  from 
the  back  side  of  the  absorber  escapes  through  the  gap  to 
the  environment.  This  limit  is  shown  in  the  dash-dotted  line 
in  Fig.  3b.  The  dashed  line  in  Fig.  3b  shows  the  effect  of  the 
geometry-dependent  radiation  losses  resulting  from  a  finite 
size  of  the  parallel  plates  on  the  performance  of  the  unicou¬ 
ple  STEG,  excluding  the  radiation  losses  from  the  thermo¬ 
electric  elements.  The  solid  line  shows  the  effect  of  both 
geometry-dependent  radiation  losses.  For  these  calcula¬ 
tions  the  electrical  contact  resistance  is  set  to  zero  and 
the  absorber  is  assumed  to  have  infinite  thermal  conductiv¬ 
ity.  The  absorber  size  is  set  to  911  mm2,  and  the  TEG 
cross-sectional  area,  Ateg,  is  varied  to  find  the  optimal  ther¬ 
mal  concentration  at  each  thermoelectric  element  length. 
The  results  show  that  excluding  the  geometry  dependence 
of  the  radiative  heat  losses  from  the  unicouple  STEG  can 
lead  to  an  overestimation  of  its  efficiency,  especially  for  a 
STEG  with  longer  thermoelectric  elements  when  the  radia¬ 
tion  contributions  of  the  elements  become  significant.  The 
contribution  of  the  radiative  heat  loss  from  the  thermoelec¬ 
tric  elements  is  very  small  for  thermoelectric  elements 
shorter  than  4  mm  which  permits  us  to  ignore  them  in 
the  radiative  heat  transfer  model  and  thus  tremendously 
reduce  the  computational  time.  Additionally,  this  observa¬ 
tion  is  the  necessary  justification  for  the  separation  of  the 
TEG  efficiency  from  the  rest  of  the  system  (Fig.  la)  as  sug¬ 
gested  in  this  paper  and  was  assumed  to  be  valid  by  several 
other  authors  (Telkes,  1954;  Rowe,  1981;  Chen,  2011). 
With  further  reduction  of  the  thermoelectric  element  length 
the  unicouple  STEG  efficiency  asymptotically  approaches 
the  efficiency  value  of  5.1%  of  a  STEG  with  infinite  parallel 
plates.  In  conclusion,  the  radiation  loss  contribution  of  the 
gap,  and  with  that  the  variation  of  the  radiative  heat  loss 
from  a  unicouple  STEG  as  a  function  of  thermoelectric  ele¬ 
ment  length,  must  be  included  in  order  to  make  accurate 
predictions  of  its  performance.  For  example,  increasing 
the  length  of  the  unicouple  TEG  from  0. 1  mm  to  3  mm  will 
reduce  the  STEG  efficiency  by  approximately  6%rel 
(Fig.  3b).  However,  assuming  that  in  a  final  STEG  design 
several  unicouple  STEGs  are  placed  close  to  each  other, 
this  effect  on  the  STEG  efficiency  is  significantly  reduced. 
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3.3.  Temperature  non-uniformity  in  the  solar  absorber 

Choosing  short  thermoelectric  elements  results  in  mini¬ 
mal  material  use  while  also  maximizing  the  potential  uni¬ 
couple  STEG  efficiency  (if  electrical  contact  resistance 
can  be  ignored)  by  minimizing  radiative  heat  losses.  How¬ 
ever,  short  thermoelectric  elements  require  large  heat  fluxes 
and  thus  a  large  thermal  concentration  to  maintain  a  large 
temperature  difference.  This  large  thermal  concentration 
can  result  in  a  large  radial  temperature  drop  within  the 
solar  absorber  from  its  outer  edge  to  the  thermoelectric  ele¬ 
ments  as  briefly  discussed  in  Section  2  (Fig.  2a). 

If  the  temperature  non-uniformity  in  the  solar  absorber 
is  large,  then  most  of  the  absorber  will  have  a  significantly 
higher  temperature  than  its  junction  with  the  thermoelec¬ 
tric  elements,  resulting  in  higher  radiative  heat  losses  than 
if  the  absorber  temperature  were  uniform  at  the  hot-junc- 
tion  temperature.  This  penalizes  the  thermal  efficiency 
and  thus  the  heat  flux  concentration.  Figs.  4  and  5  show 
the  temperature  distribution  in  the  solar  absorber  for  var¬ 
ious  thermal  concentrations.  In  Fig.  4a  the  thermal  concen¬ 
tration  is  changed  by  varying  the  absorber  size  while 
keeping  the  thermoelectric  element  cross-sectional  area 
the  same  (equivalent  radius  rteg  =  1.08  mm).  In  Fig.  4b 
the  solar  absorber  size  is  kept  constant  (equivalent  radius 
Rabs  =  17.03  mm)  while  the  cross-sectional  area  of  the  ther¬ 
moelectric  elements  is  varied.  The  temperature  drop  within 
the  absorber  in  the  latter  case  is  significantly  smaller  than 
in  the  case  with  constant  thermoelectric  cross-sectional 
area  and  changing  absorber  size.  This  is  because  the  heat 
flux  conducted  radially  is  proportional  to  the  square  of 
the  absorber  radius,  but  the  increase  of  the  radial  conduc¬ 
tion  resistance  scales  with  ln(Rabs/rteg)  (Mills,  1999).  There¬ 
fore  it  is  best  to  adjust  the  thermal  concentration  via  a 
decrease  in  cross-sectional  area  of  the  thermoelectric  ele¬ 
ments  to  keep  the  temperature  drop  within  the  absorber 
as  small  as  possible.  Using  a  significantly  larger  cross-sec¬ 
tional  area,  e.g.  with  an  equivalent  radius  of  4.8  mm  (cor¬ 
responding  to  a  TEG  module  with  20  closely-packed 
thermocouples  with  equivalent  radius  of  1.08  mm  each), 
results  in  a  detrimentally  high  temperature  drop  within 
the  absorber  in  order  to  support  the  large  flux  (Fig.  4c). 

In  addition  to  the  temperature  gradient  away  from  the 
TEG  module,  there  will  also  be  a  significant  temperature 
non-uniformity  over  the  TEG  module  (Fig.  5)  if  a  large 
number  of  thermoelectric  elements  are  closely  packed  and 
if  due  to  cost  considerations  the  copper  substrate  thickness 
t  is  chosen  to  be  small  (e.g.  0.2  mm).  However,  comparing 
Fig.  4c  and  Fig.  5a,  the  temperature  drop  in  the  TEG  junc¬ 
tion  region  is  rj5  times  smaller  than  in  the  absorber  region 
away  from  the  TEG.  This  can  be  explained  from  the  rapid 
decrease  of  the  radial  heat  flux  within  the  solar  absorber 
due  to  the  heat  removed  by  the  TEG.  However,  in  not  care¬ 
fully  designed  STEGs  this  non-uniformity  of  the  absorber 
temperature  forces  a  large  portion  of  the  TEG  to  operate 
at  off-optimal  hot-junction  temperatures  which  leads  to  a 
reduction  in  the  STEG  efficiency.  The  above  concerns 


Radial  position  on  absorber  (mm) 


Fig.  4.  Non-uniformity  of  the  solar  absorber  temperature  obtained  from  a 
radial  finite-difference  scheme.  Solar  absorber  substrate  is  copper 
Ik  =  380  W/mK)  of  thickness  t  =  0.2  mm.  Simulation  results  for  various 
thermal  concentrations,  Cth  =  with  (a)  fixed  thermoelectric  element 
radius  (rteg),  (b)  fixed  absorber  radius  (Rabs),  and  (c)  large  rteg  (fixed) 
equivalent  to  20  closely-packed  thermoelectric  unicouples. 


suggest  that  conventional  TEG  modules  with  a  large  num¬ 
ber  of  closely-spaced  thermoelectric  elements  are  less  suit¬ 
able  for  solar  applications  with  large  thermal 
concentration.  The  best  performance  will  be  achieved  for 
unicouple  STEGs  that  individually  consist  of  solar  absorb¬ 
ers  mounted  to  small  thermoelectric  unicouples  with  small 
cross-sectional  areas  (Fig.  la).  For  all  simulated  cases  up  to 
a  thermal  concentration  of  1300  the  expected  radial  tem¬ 
perature  increase  within  the  solar  absorber  away  from  the 
unicouple  TEG  is  less  than  4  °C.  The  temperature  drop 
within  the  unicouple  TEG  region  can  be  estimated  to  be 
significantly  smaller  than  1  °C.  Consequently,  it  is  reason¬ 
able  to  assume  a  uniform  absorber  temperature  equal  to 
the  hot-junction  temperature  of  the  thermoelectric  unicou¬ 
ple.  Fig.  3c  shows  the  effect  of  the  temperature  non-unifor¬ 
mity  of  the  solar  absorber  on  the  maximum  unicouple 
STEG  efficiency.  The  two  curves  show  the  difference 
between  optimizing  the  thermal  concentration  by  changing 
Ateg  versus  by  changing  Aahs  and  fixing  Ateg.  The  effect  on 
the  efficiency  as  a  function  of  chosen  element  length,  L,  is 
significantly  larger  for  the  case  with  a  constant  cross-sec¬ 
tional  area  of  the  thermoelectric  elements.  This  confirms 
that  if  the  element  length  is  to  be  reduced  in  order  to  min¬ 
imize  material  cost  without  significantly  affecting  the 
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Fig.  5.  Temperature  distribution  within  the  solar  absorber  in  the  TEG  section  with  radius  rteg  of  4.8  mm,  (a)  for  absorber  thickness  t  =  0.2mm  and 
different  thermal  concentration  C,*  and  (b)  for  Cn,  =  450  and  different  solar  absorber  thicknesses. 


uniouple  STEG  performance,  a  smaller  cross-sectional  area 
of  the  thermoelectric  elements  must  be  chosen  to  increase 
Cth. 


3.4.  Simplified  STEG  model 


If  the  contact  resistance,  the  temperature  non-unifor¬ 
mity  of  the  absorber,  the  radiation  losses  from  the  leg, 
and  the  geometry  dependence  of  the  view  factors  can  be 
ignored,  then  the  STEG  model  becomes  much  simpler  for 
three  reasons.  First,  the  absorber  efficiency  is  independent 
of  geometry.  Second,  the  discretized  TEG  submodel  can 
be  solved  independently  of  the  absorber  submodel  as  long 
as  the  temperature  and  heat  flux  at  the  junction  are 
matched.  Third,  once  a  maximum  efficiency  has  been  found 
for  given  ambient  conditions,  the  system  geometry  can  be 
scaled  without  affecting  the  efficiency  as  long  as  four 
parameters  -  Copt,  CthL,  the  length-weighted  current  iL, 
and  the  ratio  of  thermoelectric  leg  cross-sectional  areas 
AjAp  -  are  held  constant.  To  show  this  independence  of 
individual  unicouple  STEG  geometries,  the  efficiency  with 
the  above  simplifying  assumptions  can  be  expressed  as: 


=  gabs  -  grades  ~  qb.cpp 

Copfisol 

1  - 


i pL^  i„l)^abs  4rad,ss  4b,app)(CthL) ^ 

m  ' 


\1  H-  Aabs,n/ Aabstf, 

D"=(1+I(i+V/V 


(11) 

(12) 

(13) 


The  term  in  front  of  the  parentheses  in  Eq.  (11)  is  the 
optothermal  efficiency  of  the  unicouple  STEG  which  is 
often  split  into  an  optical  and  thermal  concentrating  collec¬ 
tor  efficiency  (Rowe,  1981;  Duffie  and  Beckman,  2006; 
Kraemer  et  al.,  2011;  Chen,  2011).  It  is  the  fraction  of 
the  intercepted  solar  radiation  power  which  is  delivered 
to  the  TEG.  The  heat  loss  is  determined  by  the  net  radia¬ 
tive  heat  transfer,  qmd,ss,  from  the  selective  surface  to  the 


blackbody  surroundings  (Section  2)  and  by  the  net  radia¬ 
tive  heat  loss  from  the  back  side  of  the  solar  absorber, 
kb,aPp  =  Gffd{T\s  —  TAx)  with  e,.ff  as  the  effective  emittance 
which  is  assumed  to  be  independent  of  the  thermoelectric 
element  length.  As  discussed  in  Section  3.2,  this  assumption 
introduces  the  largest  error  to  the  model,  however,  the 
inaccuracy  is  significantly  reduced  if  a  STEG  with  a  large 
number  of  closely-spaced  unicouple  STEGs  is  considered. 

The  second  term  in  the  STEG  efficiency  is  the  TEG  effi¬ 
ciency  derived  from  thermoelectric  potentials,  0  =  1/ 
u  +  ST,  with  the  relative  current  density  u  =  i/(—kS7T) 
defined  as  the  ratio  of  the  electrical  current  density  to  the 
conducted  heat  flux  (Snyder  and  Ursell,  2003).  The  model 
can  be  further  simplified  by  using  the  TEG  efficiency  equa¬ 
tion  based  on  the  averaged  ( ZT)avg  (Ioffe,  1957)  if  the  ther¬ 
moelectric  leg  is  self-compatible  (Snyder  and  Ursell,  2003). 
However,  the  following  results  include  the  temperature 
dependence  of  the  thermoelectric  material  properties  via 
the  previously  described  numerical  model  (Section  2).  Chen 
(2011)  showed  mathematically  that  CthL,  the  product  of 
the  thermal  concentration  and  the  length  of  the  thermo¬ 
electric  elements,  is  the  optimization  parameter  of  a  STEG 
if  the  thermoelectric  material  properties  are  constant 
amongst  other  assumptions.  However,  assuming  constant 
material  properties  for  the  thermoelectric  elements  and 
the  solar  absorber  can  lead  to  unrealistic  results,  signifi¬ 
cantly  over-predicting  the  conversion  efficiency  of  a  STEG 
(Chen,  1996).  Eq.  (11),  which  holds  for  temperature- 
dependent  material  properties,  shows  that  the  efficiency  is 
not  dependent  on  individual  dimensions  of  the  unicouple 
STEG  but  is  determined  by  the  single  geometric 
optimization  parameter,  CthL,  as  long  as  the  optical  effi¬ 
ciency  and  optical  concentration  does  not  change.  If  the 
thermoelectric  elements  exhibit  a  difference  in  their 
properties  of  n-  and  p-type  materials,  the  ratios  of  the  ther¬ 
moelectric  element  cross-sectional  areas,  Ap  and  An,  and  of 
the  corresponding  absorber  sections,  Aabsp  and  Aabsn,  that 
set  the  heat  input  into  the  individual  n-  and  p-type  unicou¬ 
ple  STEG  sections  (Eqs.  (12)  and  (13))  deviate  from  1  and 
are  uniquely  defined  for  the  point  of  maximum  efficiency  at 
a  specific  CthL.  Furthermore,  the  optimum  length-weighted 
electrical  current  densities,  ipL  and  inL,  are  uniquely 
defined  through  the  same  CthL.  The  relative  current 
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densities  and  with  that  the  thermoelectric  potentials  at  the 
cold  junction  are  determined  by  the  length- weighted  electri¬ 
cal  current  densities  and  therefore  are  indirectly  only 
dependent  on  Ctl,L.  In  summary,  there  is  an  optimal  CthL 
that  yields  maximum  unicouple  STEG  efficiency  with  cor¬ 
responding  area  ratios,  length-weighted  electrical  current 
densities  and  thermoelectric  potentials  if  the  STEG  is 
designed  in  such  a  way  that  the  necessary  assumptions  dis¬ 
cussed  in  this  paper  hold.  This  permits  the  length  and 
cross-sectional  area  of  the  TEG  as  well  as  the  area  of  the 
solar  absorber  to  be  chosen  arbitrarily  within  the  estab¬ 
lished  limits  without  affecting  the  STEG  performance  as 
long  as  the  optimized  CthL,  length-weighted  electrical  cur¬ 
rent  densities  and  the  area  ratios  (Eqs.  (12)  and  (13))  yield¬ 
ing  maximum  efficiency  are  unaltered.  Consequently,  it  is 
possible  to  reduce  the  volume  of  the  thermoelectric  ele¬ 
ments  and  with  that  the  thermoelectric  material  cost  for  a 
fixed  size  of  the  solar  absorber. 

4.  Applications 

4. 1.  Rooftop  STEG  performance 

So  far  we  have  discussed  the  dependence  of  the  optimal 
geometric  parameter,  CthL,  on  the  individual  dimensions  of 
the  STEG  for  fixed  cold-junction  temperature  and  incident 
solar  radiation.  In  reality  those  parameters  are  not  con¬ 
stant  over  the  course  of  a  day  or  a  year  and  are  location 
specific.  The  STEG  performance  (Kraemer  et  al.,  2011)  is 
affected  by  the  incident  solar  radiation  flux  and  cold-junc¬ 
tion  temperature  and  so  will  be  the  optimal  geometric 
parameter,  CthL. 

A  STEG  is  optimized  in  this  section  for  earth-based  elec¬ 
tricity  production  and  its  performance  is  predicted  using 
the  developed  model  neglecting  electrical  contact  resis¬ 
tances,  assuming  a  uniform  solar  absorber  temperature 
equal  to  the  hot-junction  temperature,  and  ignoring  the 
thermoelectric  elements  in  the  radiative  heat  balance.  The 
solar  absorber  is  chosen  to  be  50  mm  wide  and  of  infinite 
length  (for  view  factor  calculations);  Each  thermoelectric 
element  is  1.6  mm  long,  and  the  total  cross-sectional  area 
of  one  unicouple  TEG  is  3.645  mm2. 

The  thermoelectric  efficiency  of  the  STEG  is  strongly 
affected  by  the  amount  of  current  that  flows  through  the 
electrical  circuit.  Fig.  6a  shows  the  calculated  output  volt¬ 
age  and  efficiency  as  a  function  of  the  current  (up  to  the 
short-circuit  current)  for  a  STEG  with  CthL  =  0.4  m  with 
a  cold-junction  temperature  of  25  °C.  The  I-V curve  clearly 
shows  a  series-resistance  limited  behavior  of  the  STEG  cell 
which  is  not  due  to  resistances  at  the  electrical  contacts  or 
in  the  electrical  circuit  but  due  to  the  high  resistivity  of 
Bi2Te3  itself.  Additionally,  the  absorber  temperature  and 
therefore  the  Seebeck  voltage  decreases  with  increasing  cur¬ 
rent  which  also  affects  the  voltage  output  of  the  STEG. 

Fig.  6b  shows  the  peak  (current-optimized)  thermoelec¬ 
tric  efficiency  as  a  function  of  the  geometric  parameter  CthL 
for  various  cold-junction  temperatures.  Although  the 


optimal  geometric  parameter  is  a  function  of  the  cold-junc¬ 
tion  temperature,  the  peak  efficiency  curves  at  a  given  cold- 
junction  temperature  are  rather  flat.  At  a  cold-junction 
temperature  of  25  °C,  the  maximum  thermoelectric  effi¬ 
ciency  is  5%  at  CthL  =  0.4  m;  at  a  cold-junction  tempera¬ 
ture  of  100  °C,  the  maximum  thermoelectric  efficiency  is 
2.8%  at  CthL  =  0.35  m.  The  electrical  power  output  of  a 
unicouple  STEG  with  a  cold-junction  temperature  at 
25  °C  is  50  mW.  For  the  chosen  thermoelectric  elements’ 
dimensions  and  a  current  price  of  approximately  $300/kg 
(Bi2Te3)  the  material  cost  can  be  estimated  as  $0.27/Wp. 
Further  shrinking  of  the  unicouple  TEG  to  a  length  of 
0.5  mm  and  to  the  corresponding  cross-sectional  area  of 
~1.14mm“  to  keep  the  optimal  CthL  for  the  same  maxi¬ 
mum  STEG  power  output  reduces  the  material  cost  to 
~$0.03/Wp.  Even  the  smaller  thermoelectric  elements  are 
of  a  size  that  they  can  be  produced  with  inexpensive  bulk 
manufacturing  processes.  Therefore  the  main  cost  of  the 
unicouple  TEGs  will  be  the  thermoelectric  material. 

The  variation  of  the  incident  solar  radiation  power  not 
only  affects  the  efficiency  of  the  STEG,  but  also  the  optimal 
value  of  the  geometric  parameter  CthL  (Fig.  6c).  The  opti¬ 
mal  geometric  parameter  ( CthL)opt  decreases  (1.3  -  0.2  m) 
with  increasing  solar  insolation  (0.25  —  3  kW/m2).  How¬ 
ever,  if  the  system  is  optimized  for  highest  daily  or  all¬ 
year-round  performance  the  ( CthL)opt  will  be  slightly  differ¬ 
ent  due  to  the  daily  and  yearly  variation  of  the  incident 
solar  power.  In  order  to  optimize  a  STEG  for  maximum 
daily/all-year-round  performance  the  accumulated  electri¬ 
cal  energy  output  of  the  STEG  must  be  maximized. 
Fig.  7  shows  (a)  the  simulated  solar  intensity  distribution 
for  Boston  over  the  course  of  a  simulated  sunny  equinox 
day  (September  22nd)  (Bird  and  Riordan,  1984)  and  (b) 
the  accumulated  electrical  energy  output  per  day  of  a  cur¬ 
rent-optimized  STEG  as  a  function  of  the  CthL.  In  those 
calculations  the  absorptance  of  the  solar  absorber  is 
assumed  to  be  independent  of  the  incident  angle  of  the 
solar  radiation.  For  incident  angles  up  to  60°,  this  is  a  rea¬ 
sonable  assumption  for  solar  absorbers  based  on  multi¬ 
layer  selective  surfaces  which  are  most  commonly  used 
for  solar  thermal  applications  (Tesfamichael  and  Waeckel- 
gard,  2000).  Additionally,  possible  shading  effects  occur¬ 
ring  in  a  final  design  structure  are  neglected.  For  final 
device  performance  predictions  those  effects  must  be 
included.  Independent  of  those  details,  the  daily  electrical 
performance  of  the  STEG  only  depends  weakly  on  C,hL 
in  a  wide  range  (0.35  —  0.8  m)  around  the  maximum  of 
0.3  kWh/m2/day.  Consequently,  within  that  range  the  CthL 
can  be  chosen  for  highest  overall  system  efficiency  if  the 
STEG  is  integrated  for  example  in  a  cogeneration  or  hybrid 
system  (Kraemer  et  al.,  2008)  without  significantly  affecting 
the  STEG  performance. 

4.2.  Rooftop  electric  power  and  hot  water 

In  solar  thermoelectric  systems,  approximately  50-80% 
of  the  intercepted  solar  heat  is  released  at  the  cold  junction 
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Fig.  6.  STEG  performance  characteristics,  (a)  Voltage  and  power  output  as  a  function  of  current  for  an  optimized  unicouple  STEG  (C,i,L  =  0.4  m)  under 
1  kW/m2  illumination,  cold-junction  temperature  of  25  °C,  chosen  thermoelectric  element  dimensions  (W  x  D  x  L)  of  1 .35  mm  x  1 .35  mm  x  1 .6  mm,  and 
1  m2  area;  (b)  STEG  efficiency  as  a  function  of  the  geometric  parameter  CthL  for  different  cold-junction  temperatures;  (c)  STEG  efficiency  as  a  function  of 
incident  solar  intensity  for  different  ClhL:;  and  (d)  excess  heat  released  at  the  cold  side  for  various  cold-junction  temperatures  as  a  function  of  the  STEG 
current. 


Time  C^L  (m) 


Fig.  7.  STEG  optimization  over  the  course  of  a  day.  (a)  Simulated  solar  intensity  distribution  for  Boston  over  the  course  of  the  day  of  September  22nd; 
and  (b)  the  accumulated  electrical  energy  output  per  day  of  a  current-optimized  STEG  as  a  function  of  the  CthL. 


of  the  device  (Fig.  6d).  This  waste  heat  must  be  removed  in 
order  to  maintain  the  cold  junction  at  a  given  temperature 
for  the  highest  STEG  performance.  One  way  to  remove  this 
excess  heat  is  first  to  spread  it  out  on  the  cold  side  using  a 
metallic  heat  spreader,  and  then  to  reject  the  heat  to  the 
environment  via  natural  convection,  similar  to  what  is 
done  for  the  heat  management  of  PV  cells.  In  certain  appli¬ 
cations  such  as  cogeneration,  this  waste  heat  is  actually  the 
input  to  a  secondary  system  such  as  a  domestic  hot  water 
loop.  Small  deviations  from  the  peak  TEG  operating  point 
can  have  large  (positive  or  negative)  effects  on  the  quantity 
of  heat  delivered  to  the  heat  sink.  As  a  result,  a  STEG 
cogeneration  system  can  be  designed  to  favor  the  produc¬ 
tion  of  either  electrical  power  or  waste  heat,  depending 
on  the  demands  of  the  application.  This  can  be  accom¬ 
plished  by  choosing  the  appropriate  geometric  parameter 


CthL  (Fig.  6b  and  d).  If  it  is  desired  to  optimize  the  system 
for  maximum  electrical  power  then  there  is  one  specific 
optimal  (ClhL)opt  for  a  specific  cold-junction  temperature 
and  incident  solar  flux  as  discussed  in  previous  sections. 
Conversely,  choosing  a  smaller  CthL  results  in  more  heat 
transported  to  the  hot  water  loop  because  the  absorber 
stays  at  a  lower  temperature  and  thus  the  radiative  heat 
losses  of  the  system  are  smaller.  Fortunately,  as  mentioned 
in  the  previous  section,  the  daily  total  STEG  electrical 
energy  output  is  a  weaker  function  of  CthL,  so  it  is  possible 
to  deviate  from  the  geometric  parameter  optimized  for 
maximum  electrical  performance  with  only  a  negligible 
effect  on  the  daily  electrical  performance. 

For  a  system  with  fixed  CthL,  it  is  even  possible  to  adjust 
the  balance  between  electrical  power  and  waste  heat  solely 
by  adjusting  the  electrical  current  of  the  circuit.  As  an 
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STEG  cold-junction  temperature  (°C) 


Fig.  8.  STEG  operation  diagram:  excess  heat  available  for  cogeneration  at 
the  cold  junction  of  a  STEG  (optimized  for  1  sun  illumination,  and  cold- 
junction  temperature  of  25  °C)  of  1  m2  area  under  1  kW/m2  illumination 
as  a  function  of  the  cold-junction  temperature.  Also  shown  are  the  STEG 
efficiency  contour  lines  (dashed),  STEG  cell  iso-current  lines  (solid),  the 
iso-overall-heat-transfer-coefficient  ( U)  lines  (dash-dotted),  and  the  iso¬ 
therm  lines  of  the  solar  absorber  (dotted). 

example,  we  consider  a  system  of  CthL  =  0.4  m.  Because 
the  available  excess  heat  and  cold-junction  temperature 
are  both  affected  by  the  current,  it  is  useful  to  plot  the  heat 
available  for  cogeneration  per  square  meter  of  absorber 
area  as  a  function  of  cold-junction  temperature  at  various 
currents  (Fig.  8).  In  those  calculations  the  possible  radia¬ 
tive  heat  losses  from  the  back  side  of  the  copper  heat  sink 
to  the  glass  enclosure  (Fig.  la)  are  taken  into  account  and 
reduce  the  available  heat  for  cogeneration.  Superimposed 
on  Fig.  8  are  two  sets  of  lines:  red2  dotted  lines  are  contour 
lines  representing  the  corresponding  STEG  efficiency;  blue 
dotted  lines  are  the  corresponding  absorber  temperature. 
This  STEG  operation  diagram  can  be  used  to  determine 
the  conditions  where  the  STEG  can  operate,  because  the 
performance  of  any  heat  sink  can  be  characterized  by  its 
relationship  between  the  heat  sink  temperature  and  the 
rejected  heat.  The  constant  of  proportionality  between 
the  transferred  heat  and  the  temperature  difference  is  the 
overall  heat  transfer  coefficient,  U  in  units  W/K  of  the  heat 
sink.  In  Fig.  8,  dash-dotted  lines  of  constant  U  are  plotted 
in  green  for  the  example  of  a  heat  sink  comprising  a  25  °C 
fluid  passing  over  the  cold  junction.  Changing  the  electrical 
current  allows  the  system  to  operate  at  different  points 
along  this  heat  transfer  coefficient  characteristic  curve, 
affecting  the  cold-junction  temperature,  the  amount  of 


2  For  interpretation  of  color  in  Figs.  1-3  and  6-8,  the  reader  is  referred 
to  the  web  version  of  this  article. 


rejected  heat,  and  the  amount  of  electrical  power  gener¬ 
ated.  For  example,  if  it  is  desirable  to  generate  more  waste 
heat  in  the  morning  and  more  electricity  in  the  afternoon,  it 
is  possible  to  run  the  system  in  an  “overdrive”  mode  (with 
super-optimal  current)  in  the  morning,  and  then  run  the 
system  at  optimal  current  in  the  afternoon.  It  is  interesting 
to  notice  that  there  will  always  be  two  operational  current 
points  with  the  same  STEG  efficiency  but  with  different 
cold-junction  heat  rejection  rates  if  the  STEG  is  operated 
at  off-optimal  conditions.  For  example  for  a  heat  sink  with 
overall  heat  transfer  coefficient  of  25  W/K  the  line  of  con¬ 
stant  U  intersects  the  4%  efficiency  contour  line  twice.  One 
intersection  corresponds  to  a  current  of  approximately 
0.8  A  and  the  other  one  to  1.25  A.  At  the  higher  STEG  cur¬ 
rent  the  effective  conductance  of  the  TEG  increases  due  to 
thermal  transport  via  electrons  and  holes,  which  results  in  a 
lower  absorber  temperature  (198  °C).  This  decreases  the 
radiation  losses  but  also  decreases  the  TEG  conversion  effi¬ 
ciency.  These  two  effects  result  in  less  electricity  generation, 
and  more  heat  transported  to  the  cold  side  (640  W).  In 
addition,  the  cold-junction  temperature  increases  slightly 
because  of  the  assumed  constant  U-value.  At  the  lower 
STEG  current  the  effective  conductance  of  the  TEG 
decreases,  resulting  in  a  slightly  lower  cold-junction  tem¬ 
perature,  a  higher  absorber  temperature  (211  °C),  and 
higher  radiation  losses,  but  a  better  TEG  efficiency.  The 
net  result  is  the  same  4%  STEG  efficiency,  yet  in  this  case 
less  heat  is  transported  to  the  cold  side  (605  W). 

5.  Conclusion 

We  developed  a  model  and  optimization  methodology 
for  earth-based  STEG  designs  using  a  glass  vacuum  enclo¬ 
sure  and  a  large  flat-plate  thermal  concentration.  A  few 
assumptions  are  made  and  validated  to  present  a  simplified 
yet  useful  and  realistic  model  for  researchers  and  engineers 
who  want  to  design  STEGs  or  want  to  incorporate  them  in 
cogeneration  systems.  We  describe,  discuss  and  justify  the 
necessary  assumptions  that  have  to  hold  for  the  STEG 
design  to  be  independent  of  individual  dimensions.  The 
only  geometric  parameter  which  has  to  be  optimized  for 
highest  STEG  efficiency  is  CthL,  even  if  the  temperature 
dependence  of  the  thermoelectric  materials  and  the  solar 
absorber  is  included.  This  is  a  very  useful  result  because 
it  shows  that  the  performance  of  a  STEG  is  independent 
of  the  thermoelectric  material  volume,  and  the  cost  for 
the  thermoelectric  material  can  be  negligibly  small  com¬ 
pared  to  the  total  system  cost.  The  numerical  simulations 
revealed  that  the  optimal  ( CthL)opt  is  weakly  dependent 
on  the  cold-junction  temperature,  but  shows  a  stronger 
dependence  on  the  incident  solar  radiation.  If  the  STEG 
efficiency  is  averaged  over  the  course  of  a  day  only  small 
variations  in  the  efficiency  are  observed  in  a  wide  range 
of  CthL  around  the  maximum.  Consequently,  if  a  STEG 
is  integrated  in  a  cogeneration  system  the  CthL  can  be  cho¬ 
sen  within  that  range  to  optimize  the  overall  efficiency 
without  significantly  affecting  the  STEG  performance. 
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For  the  matter  of  heat  removal  at  the  cold  junction  of  the 
STEG  we  discussed  the  influence  of  the  current  on  the 
released  excess  heat  and  the  STEG  efficiency  and  intro¬ 
duced  a  STEG  operation  diagram.  This  diagram  illustrates 
the  possible  region  of  operation  for  the  optimized  STEG 
with  a  specific  heat  sink  fluid  temperature  and  solar  insola¬ 
tion.  The  current  is  an  additional  parameter  to  actively 
optimize  the  operation  of  the  cogeneration  system  and 
therefore  the  model  underlying  the  STEG  operation  dia¬ 
gram  can  be  used  in  controllers  of  those  systems.  The  pre¬ 
dicted  efficiency  for  the  Bi2Te3  based  STEG  considered  in 
this  paper  is  5%  at  AM1.5G  conditions  and  with  a  cold- 
junction  temperature  of  25  °C.  With  this  efficiency  and 
the  concept  of  large  thermal  concentration  the  thermoelec¬ 
tric  material  cost  can  easily  be  reduced  to  below  $0.05/Wp 
and  become  negligible  compared  to  the  overall  system  cost. 
Consequently,  integrating  STEGs  in  solar  hot  water  vac¬ 
uum  tube  systems  for  cogeneration  application  shows 
potential  to  provide  cost-competitive  domestic  hot  water 
and  electricity.  However,  more  research  and  development 
is  needed  to  explore  this  potential. 
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